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ABSTRACT 



Aims. We report a serendipitous detection of rapid, large amplitude flux density variations in the highly core-dominated, flat-spectrum 
radio quasar 1 156+295 during an observing session at the Very Long Baseline Array (VLBA). 

Methods. The source was observed as a part of the MOJAVE survey programme with the VLBA at 15 GHz on February 5, 2007. 
Large amplitude variability in the correlated flux density, unexplainable in terms of the source structure, was first discovered while 
processing the data, and later confirmed by calibrating the antenna gains using 24 other sources observed in the experiment. 
Results. The source shows variations in the correlated flux density as high as 40% on a timescale of only 2.7 hours. This places 
1 156+295 between the classical IDV sources and the so-called intra-hour variables. The observed variability timescale and the mod- 
ulation index of 13% are consistent with interstellar scintillation by a nearby, highly turbulent scattering screen. The large modulation 
index at 15 GHz implies a scattering measure that is atypically high for a high galactic latitude source such as 1 156+295. 

Key words. Galaxies: active - Galaxies: jets - quasars: individual: 1 156+295 - Techniques: interferometric 



1. Introduction 

Intraday variability (IDV) of compact, radio-loud AGN at cm- 
wav elengths was first discovered in the mid-1980s (■Witzel et alJ 
11986: .Heeschen et al.l IT987i) and both source-intrinsic and ex- 
trinsic mechanisms for IDV have been vigo rously studied 
(for a review, see e.g. IWagner& Witzeilll995h . Detection of 
time delays in the variability pattern arrival times between 
widely separated telescopes, as well as observed annual mod- 
_ ulation of the variability timescale, have shown conclusively 
, that interstellar scintillation (ISS) is the cause of intra-hour 
' flux density variations observed in the three most extreme 
; IDV sources PKS 040 5-385 (Jauncev et al] l2000l) . J1819+3845 
" (Dennett-Thorpe & de Bruvn 20021 120031) . and PKS 1257-326 
(.BignaU et al. 2003. .2006,) ! Evidence for the ISS or igin of IDV 
was found also in the cases of 0917+624 (Ricket t etal J 120011; 
iJauncev & Macquart' '2001': iFuhrmann et alj |2002|). PKS 1519- 
273 (Jauncev et al. 2003), and Jl 128+5925 dGabanvi et"an 
|2007|) . This provided the possibility to study both small-scale 
spatial fluctuations in the interstellar medium and the struc- 
ture of compact radio sources on the microarcsecond scale. 
Unfortunately, the extreme IDV sources, which are the most suit- 
able for these studies, are rare. By "extreme", we mean sources 
showing variability on timescales of a few hours or less and with 
an rms amplitude of modulation of over 10%. In this Letter, we 
report the serendipitous discovery of very large amplitude IDV 
in quasar 1156+295 with a timescale of variations shorter than 
3 h. This discovery is also unusual because it was achieved by a 
VLBI experiment. 

1156+295 (4C + 29.45) is an optically violently variable 
quasar at z = 0.729 dBurbid ge 1968), which has Galactic coor- 
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Fig.l. Naturally weighted 15 GHz VLBA image of 1156+295 
observed on February 5, 2007. The map peak brightness is 
1.46 Jy beam"'. The contours begin at 0.7mJybeam"' and in- 
crease in steps of 2. Since the source flux density varied signifi- 
cantly during the observation, an amplitude self-calibration had 
to be used early in the imaging process. This may affect the im- 
age quality. 



dinates / - 199.4°, b - +78.4°. The source is strongly variable 
throughout the electromagnetic spectrum from radio to y-rays. 
At radio frequencies, 1156+295 shows significant long-term 
variability on a timescale of months (e.g. iKovalev et al.ll2002h . 
Wills et al. ( 1983) reported a range of at least 5 magnitudes in the 
optical brightness, and in gamma-rays ISreekumar et alj (Il996h 
detected an order of magnitude variability using the EGRET 



2 



T. Savolainen and Y. Y. Kovalev: Large-amplitude, intraday variability in 1156+295 




8 J 



in o 

IF2 + 

IF3 o 

IF4 " 




in □ 

ira + 

IF3 o 

IF4 A 



la" le"* so' i''o'' o'd' 



♦ 

! * 



* ♦ 



13° 16° 20° ro° 





MK-SC 















Time 



Time 



Time 



Fig. 2. Top: Amplitude self-calibration solutions of three example antennas (Brewster, Mauna Kea, and Pie Town) for the VLBA 
experiment on February 5, 2007. There is one solution per scan, and IFs are kept separate. Low elevation scans (< 15°) are excluded 
from the figure. Red symbols correspond to scans of 1156+295 and blue symbols correspond to all other sources. The obvious 
discrepancy between the solutions for 1156+295 and the other sources is present at all ten antennas. Bottom: Examples of the 
correlated flux density curves of 1156+295 at three baselines after calibration of antenna gains by amplitude correction factors 
interpolated from the nearby scans of other sources. Time is in UT. 



instrument onboard the Compton Gamma Ray Observatory. 
During high radio brightness states, 1156+295 appears to be 
highly core-dominate d in VLBI images (95% at some epochs 
jKovalev et al.l l2005l) : see also Fig. [T]i. Superluminal motion 
at s peeds ranging from 3.5/i ~ 'c to 14.1/i~' c were reported by 
e.g. iPiner & K inghamI ([T997h . iJorstad et al.l (1200 ih . iHong et al.l 
(|2004|), andlKellermannetaD (l2004l) . 

In the optical bands, 1156+295 shows rapid variations 
on timescales of between less than 30 minutes and 3 days 
dWifls et al.lll"983t iRaiteri et al.lll998b . The first I DV detection 
of 1 15 6+295 at radio frequencies was reported by iLovell et"an 
(I2003h . who observed the source with the VLA at 5 GHz in 
January 2002. Their study measured rms flux density variations 
of 167 mJy during 3 nights of observations and a variability 
timescale of ~ 24 h. The re lative amplitude of variations ob- 
served bv lLovell et al.l ( l2003l) was far smaller than reported here. 

Throughout this paper, we use the following cosmological 
parameters: Hq = 73 km s"' Mpc"', - 0.27, and Qa = 0.73. 
The angular scahng conversion for z - 0.729 is 7.07 pc mas ' . 

2. Observations and analysis 

Quasar 1 156+295 and 24 other compact radio sources were ob- 
served with the NRAO's Very Long Bas eline Array (VLBA) at 
15 GHz as a part of the MOJAVE project (iLister & Homanl2005h 
on Februai-y 5, 20070. The entire observing session had a dura- 
tion of 24 hours, which included 9 scans of 1 156+295, each last- 
ing 4.7 minutes. The signal was recorded in dual-polarisation 
mode: 4 IFs (each with 8 MHz bandwidth) per circular po- 
larisation and 1-bit sampling. The data from the experiment 
were calibrated fo llowing st andard procedures of VLBI data re- 
duction (see e.g. .Lister & Homan 2 005) . After a priori ampli- 

' We note that 1 156+295 was observed several times in the MOJAVE 
survey and its IDV behaviour during the rest of the MOJAVE epochs 
will be analysed in Kuchibhotla et al. (in prep.). 



tude calibration and fringe fitting, the correlated flux density of 
1156+295 showed significant, correlated, temporal variability 
from scan to scan on every baseline: between 7h and 12h UT, 
it first dropped by about 0.6 Jy and then increased again. The 
dip can be seen on every baseline, which excludes the source 
structure as the cause of the variations. The system temperature 
measurements used in the a priori amplitude calibration do not 
show any jumps that could explain the dip. 

To ascertain whether the variability is genuine or due to am- 
pUtude calibration problems, we first imaged and self-calibrated 
the (m, v) data of 1 156+295 in a normal manner and then com- 
pared the resulting antenna gain correction factors with those de- 
rived from the self-calibration of the other 24 sources observed 
in the same session. Amplitude self-calibration with a 5-minute 
solution interval was indeed able to remove the variability of 
1 156+295, but it produced gain amplitude correction factors that 
were -for every antenna - significantly offset from the gain cor- 
rections determined from the 24 other sources in the experiment 
(see examples in Fig.|2]i. This can be explained if the variability 
is genuine but does not alter the shape of the source's visibil- 
ity function significantly during the experiment, e.g. if the vari- 
ability occurs in one dominant compact component. Since the 
antenna gain is not source-dependent, calibration errors cannot 
account for the above-described behaviour. 

Although improbable, there is a possibility that the observed 
variations could be due to a source-specific correlator problem. 
In such a case, however, one would expect to see discrepant solu- 
tions from the fringe fitting. We carefully checked the fringe-fit 
solutions and could not find any anomalies. Therefore, an in- 
strumental effect that originates in the correlator seems highly 
unlikely. Taken together, we conclude from the above arguments 
that the variations in the correlated flux density of 1 156+295 are 
due to genuine intraday variability observed in the source. 

To study rigorously the observed IDV event in 1156+295, 
the residual errors in the antenna gains, which remain after the a 
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Fig. 3. Calibrated coiTelated Stokes I flux density of 1156+295 
(in Jy) as a function of (u, v) radius (in MA) for each individual 
scan of the VLBI experiment on February 5, 2007. The data has 
been averaged over the IFs and for 30 s in time. The different 
scans are shown in diff'erent colours with the scan start time in 
UT indicated in the legend. 
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Fig. 4. Integrated flux density curve of 1156+295, obtained 
by averaging the correlated flux density at projected baselines 
shorter than 100 M/1. The error bars represent the standard de- 
viation of the correlated flux density (at baselines < lOOM/l) in 
each scan. Colour coding of the scans is the same as in Fig. [3] 



priori amplitude calibration, need to be corrected. Since a large 
number of sources were observed during the VLBA session with 
their scans interleaved, it was possible to estimate the antenna 
gain correction factors for the scans of 1156+295 by linearly 
interpolating the amplitude self-calibration solutions from the 
scans of the other sources that were observed within one hour 
of each scan of 1 156+295. Before interpolation, the scans at el- 
evations below 15° were removed. The bottom panels in Fig.|2] 
show examples of correlated flux density curves on three base- 
lines that were calibrated in this way. The variations in the corre- 
lated flux density are very similar at short (LA-PT), intermediate 
(BR-OV), and long (MK-SC) baselines. Again, this is expected, 
if the variability occurs in a single dominant compact compo- 
nent. In Fig. [3] we plot the calibrated correlated flux density of 
1 156+295 as a function of (u, v) radius for each individual scan. 
The figure shows that the shape of the visibility function, "Viru^v), 
changes little from scan to scan, although the correlated flux den- 
sity changes by 0.6 Jy. To analyse possible small changes in the 
shape of 'y(r„_,,), we constructed the mean visibility function of 
the observation, ('V(r„ i )), by averaging the normalised visibil- 
ities of all individual scans. The correlated flux density of each 
scan was first binned into 50 MA wide bins in projected (u, v) 
spacing and then normalised to 1 .0 at - 25 MA. We sub- 
tracted (T^(r„,i,)) from the binned ^(r„ ,,) of each individual scan 
and inspected the residuals. Some variability was found: T^(r„ i.) 
at 200 < r„ y < 300 M/l increased by 8% during the observation, 
which is more than the ^ 3% variability expected from the time 
variable response of the interferometer to the faint jet. However, 
the variability turned out to be solely due to the baselines to St. 
Croix, and therefore probably was due to calibration inaccura- 
cies of the antenna. Consequently, we are unable to confirm con- 
clusively the detection of time variability in T^(r„ i ). An upper 
limit to variability is 8%. 

We constructed an integrated VLBA flux density curve for 
1156+295 by averaging the correlated flux density at projected 
baselines between 6 - 100 M/l (Fig.|4]l. Since 'y{ru,v) was essen- 
tially flat in this range (Fig. [3]), the averaging provided a good 
estimate of the integrated emission coming from angular scales 
< 30 mas. Averaging also reduces the errors due to inaccurate 
amplitude calibration, because these are antenna-specific. 



The two important quantities in the analysis of IDV are the 
modulation index and the characteristic timescale of variabil- 
ity. We consider the average of the peak-to-trough and trough- 
to-peak times of the large dip to be the timescale of vari- 
ability, which is therefore fvar = 2.7 + 0.5 h. We note that 
our fvar is approximately a factor of 1.7 times longer than 
the timescale on which the autocorrelation function of inten- 
sity fluctuations reaches a fraction 1/e of its maximum value 
dJauncev & MacquartI 12001*). The modulation index m, defined 
to be the standard deviation of the source flux density divided 
by the mean source flux density, is 13 + 3% for the flux den- 
sity curve in Fig. H] The uncertai nties in fyar and m were es- 
timated following the analysis in iDennett-Thorpe & de Bruvnl 
(2003). The modulation index is much higher than the val- 
ues typically observed for IDV sourc e s at frequencies a bove 
5 GHz (Kedziora-Chu dczer et alj|200ll: iKraus et alj|2003l) . For 
1 156+295. iLovell et akf ^OOSb reported m = 5.8% at 5 GHz in 
2002. 



3. Discussion 

If we assumed a source-intrinsic origin of the observed vari- 
ability in 1156+295, it would imply, by light travel time ar- 
guments, a brightness temperature of > 2 x 10'^ K, which 
is far in excess of the inverse C ompton (IC) limit of 10'^ K 
(iKellermann & Paulinv-Totiil 11969 ). A Doppler factor higher 
than ~ 270 would be then required to avoid the IC catastrophe. 
This would cause severe problems with the energy requirements 
in the source (Begelman et al. 1994). Extremely fast jets are also 
unlikely in the light of VLBI monitoring surveys, wh ich indicate 
that t he maximum jet Lorentz factor F is 30 - 40 ( Cohen et al.l 
l2007h . The standard model of incoherent synchrotron radiation 
from relativistic electrons in the jet would therefore have signif- 
icant difficulties in explaining the observed variability, if it were 
intrinsic. For this reason, we concentrate on two possible ex- 
trinsic causes: interstellar scintillation and an extreme scattering 
event. 

Propagation effects in the ionised medium of the Milky Way 
will cause a suffi ciently compact source to scintillate (for review 
of ISS, see e.g. lRickettlll990t lGoodmanl ll997). Since the ob- 
served m in 1156+295 is clearly below 100%, it is reasonable 
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to assume that difFractive scintillation is fully quenched and to 
consider only refractive ISS. The observed fyai and m constrain 
the properties of a possible scattering screen and the size of the 
scintillating source. Because of the high m in our case, it ap- 
pears likely that our observing frequency is close to the criti- 
cal frequency Vs between strong and weak scattering regimes. 
iGoodmaiil (1 19971) provided a formula for m that can be used 
to interpolate solutions close to Vs, where the asymptotic solu- 
tions for strong and weak regimes break down. Combining his 
Eqs. 10, 12, 18, 19, and 20, and assuming a Gaussian bright- 
ness profile for the source, we calculated the screen distance, 
Dscr, and the FWHM size of the scintillating source O^^"'^ as a 

function of scattering measure SM - Jj^ Cj^dx, where Cjj is the 

strength of the electron-density fluctuation^ (Fig. |5]). The re- 
sults were parameterised by the screen velocity. We estimated 
the uncertainties in Dsa and 6^^"^ by Monte Carlo methods 
and the hatched regions in Fig. |5]show the Icr error for these 
quantities. In additi on, the results given by the interpolation for- 
mula of Goodman! (119971) were confirmed by numerically inte- 
grating the intensity covari ance function for r efractive scintil- 
lation (see Appendix A in iRickett et al.l l2b06h . In the bottom 
panel of Fig. |5] we indicate an approximate lower limit to the 
source size, g^^^'^ > H ij,as, given by the IC catastrophe limit to 0; 



and by assuming an upper lim i t of 50 for the Do ppler factor 
jLahteenmaki & Valtaoialll999l; ICohen et al.ll2007l) . As can be 
seen from Fig.|5] if the variability is due to ISS, we have either a 
rather nearby screen (Dsci- < 300 pc) with SM > 0.5 m^^o/^ pc or 
qFwhm small that it would require Doppler factor exceeding 
50. 

We can estimate g^^f^'^ independently by assuming that 
there is approximate equipartition between the energy den- 
sities of the rnagnetic field and the radiating electrons 
(IScott & Readheadlll977h . We assume the synchrotron peak fre- 
quency to be 15 GHz and the corresponding peak flux density to 
be 1.5 Jy. This is based on the assumption that the s cintillating 
component is the core of a i Blandford & Koni"gl (I1979D type jet 
and corresponds to the t - 1 surface at our observing frequency. 
The resulting equipartition size is g^^^'^ ~ 270 ■ 5"^^^^juas, 
where 6 is the Doppler factor. If we again take 6 < 50, this results 
in 6^"'™ > 55 //as, SM> 4m-^°'^pc, and D,^, < lOOpc for 
screen velocities between 10 and 50kms~'. S i nce th e Galactic 
electron distribution model by lCordes & Laziol (|2002b predicts a 
scattering measure of only 0. 1 m"^"^^ pc for the line-of-sight to- 
wards 1 156+295, our results indicate a nearby, localised region 
of highly turbulent ionised gas in that direction. 

The frequency-dependence of m is slightly su rprising: one 
would not expect m to increase from 5.8% at 5 GHz jLovell et alj 
l2003h to 13% at 15 GHz, unless the scattering is strong and 
the source is smaller than the scattering angle. In that case the 
scintillating component would contain only about 30 - 50% of 
the total flux density of the core, because otherwise the ob- 
served m would be significantly higher. Another possibility is 
that 1156+295 was in a more compact stage during our obser- 
vation than it was 5 years earlier. Simultaneous multi-frequency 
measurements of m are required to clarify this. 

It is also possible that the observed variations are not due 
to standard ISS, but in stead correspond to an isolated extreme 
scattering event (ESE: [Fiedler et al]|l987h . Unfortunately, our 
short, single-frequency observation does not provide us with 
sufficient information to distinguish between these two cases. 
However, the shape of the flux density curve resembles an ESE 
with two maxima symmetrically surrounding a deep minimum. 




Fig. 5. Parameters of the possible scattering screen. Top: The 
distance to the screen as a function of the scattering measure 
SM. Different screen velocities are shown in different colours. 
Bottom: The scintillating source size as a function of SM. The 
right edge of the panel gives the Doppler factor if g^^'^*^ is 
equal to equipartition size. The grey shaded area corresponds 
< 17yuas. In both panels, the hatched region around 



fWHM 



Note that 0, in lGoodmad ( Il997h is approximately (9f "'""/2.35. 



the solid lines gives Icr error for the plotted quantities. The +lcr 
lines in the g^^^'^ plot roughly correspond to -Icr lines in the 
Dscr plot and vice versa. 



Althoug h "classical" ESEs have durations of between weeks and 
months. ' Cimo et al] (|2002|) reported an event in 0954+658 with 
a duration of less than 2 days. The extremely short timescale 
of our event would, in the case of an ESE, imply a cloud size 
< 0.01 AU. 
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